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Abstract 

A nonrelativistic Hamiltonian with plausible spin dependent cor- 
rections is proposed for the quarkonia below their respective strong 
decay thresholds. With only six free parameters this model repro- 
duces the nine known masses of the bottomonia within about 1 MeV, 
the six known masses of the charmonia within a few MeV and the 
five known leptonic decay widths of the 3 S± states within about 20 %. 
The model is then used to predict the masses of the remaining 43 
quarkonia (some of them for the first time) and of the leptonic decay 
widths of the two 1 So(bc) states. Comparison with some other models 
is made. 
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1 Introduction 



In the present paper we study the mass spectra of heavy quarkonia below 
their strong decay thresholds, i.e. below 10558 MeV for the bb quarkonia, be- 
low 7144 MeV for the be quarkonia and below 3729 MeV for the cc quarkonia. 
For the quarkonia with unnatural parity, which because of the conservation 
laws cannot decay strongly into two pseudoscalar mesons, these thresholds 
should be a little higher. In practice this distinction is important only for 
the quarkonia 2Px(bc) and 2P]/(&c), which according to our calculation have 
masses above the minimal thresholds given above, but below their real thresh- 
old, which is 7189 MeV. As a byproduct we obtain the leptonic widths of the 
3 S\ states of the bb and of the cc quarkonia, as well as of the l So states of the 
be quarkonia. Experimentally, out of the 34 bb quarkonia expected nine have 
been observed. Here and in the following we consider a particle as observed, 
if it is listed as firmly established in the 1996 Particle Data Group Tables [p]]. 
Out of the eight expected cc quarkonia six have been observed and only the 
singlet P state and the excited t] c are still missing. The masses of both have 
been reported [Q , but they are not considered as firmly established. None of 
the be quarkonia has been observed, but candidates have been reported @ 
and discoveries are expected in the near future. 

In our previous paper || (further quoted I) we have pointed out that a 
simple nonrelativistic model can reproduce, among other things, the masses 
of the known 3 S(bb) states and of the centres of gravity of the known 3 P(bb) 
states within the experimental errors. After this paper had been published, 
we became aware of a series of papers (cf. ref. j| and papers quoted there) 
using essentially the same potential. There are many models, which give 
good fits to the masses of the bb quarkonia. The review paper || quotes 
and discusses about 30 of them. Our model, however, seems to be the only 
one so far, which reproduces these masses within the experimental errors, 
i.e. with a precision of about 0.5 MeV, which corresponds to 50 ppm. of the 
total mass, or to 0.1 per cent of the first excitation energy. This result is 
amazing, because the mean square velocity of a 6-quark in the ground state 
of the bb system is, in a system of units where the velocity of light c = 1, 
(v 2 ) ~ 0.08. Thus a simple-minded estimate of the relativistic corrections 
would give (v 2 ) 2 ~ 0.6% of the total mass. The usual interpretation is 
that the nonrelativistic Hamiltonian is an effective Hamiltonian and that the 
relativistic corrections are taken into account by a renormalization of the 
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parameters of this Hamiltonian. 

In the present paper we extend the model described in I in two ways. 
Firstly, we generalize the nonrelativistic potential so that it applies also to 
be and cc quarkonia. This requires one more parameter — the mass of the 
c-quark. We find that the predictions for the cc quarkonia agree with exper- 
iment within about 4 MeV i.e. within about 0.1 per cent of the total mass, 
or one per cent of the first excitation energy. Since the mean square velocity 
of the c-quark in the ground state of the cc system is about 0.25, i.e. the 
mean root square velocity is about half the velocity of light, this good fit is 
even more striking than that for the bb case. Then we make predictions for 
the yet undiscovered states hoping that, since our fits for the known states 
are good, the predictions will also work. One should keep in mind, however, 
that the problem of the be quarkonia is not just an extrapolation between the 
bb and cc cases. If our handling of the effects of the mass difference between 
the two constituents is faulty, the error of the predictions may be larger than 
expected. 

Secondly, we supplement our nonrelativistic Hamiltonian with the stan- 
dard spin-dependent terms. With one more parameter - the coupling con- 
stant a s (m 2 ) - we describe all the hyperfine and fine splittings, as well as the 
leptonic decay widths. The value of a s (m 2 ) is so plausible - it corresponds 
to a s {m 2 z ) = 0.115 - that one can either interpret it as a determination of 
a s (m 2 z ), or as a known quantity and not a free parameter. We choose the 
former possibility. By comparison with experiment we find that the errors of 
the calculated splittings do not exceed about 1 MeV for the bb systems and 
about 5 MeV for the cc systems. Also the leptonic decay widths agree reason- 
ably well with experiment. Again by analogy we expect that our predictions 
for the yet undiscovered states should be good. 



2 Spin averaged masses of spin triplet states 

The nonrelativistic potential used in the present paper for the colour-singlet 
QQ-systems, where Q may, but does not have to, be the charge conjugate of 
Q, is 

V(r) =mQ + m Q - 0.78891 + 0.70638^ - 0.32525-, (1) 
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where all the constants are in suitable powers of GeV and the quark masses 
are 



m b = 4.8030 GeV, 



m c = 1.3959 GeV 



(2) 



For the be quarkonia we use the reduced mass 



1.0816 GeV. 



(3) 



Hbc 



nib + m c 



The parameters are given with so many digits only in order to assist the 
reader, who would like to check our calculations. We do not claim to have 
established the quark masses with such accuracies. Similar remarks apply 
to coupling constants etc. given further. For the bb system this potential 
reduces to the potential given in I. 

According to our interpretation the eigenvalues of the nonrelativistic 
Hamiltonian with this potential should be interpreted as the masses of the 
centres of gravity of the spin triplets. Another popular interpretation is that 
the nonrelativistic Hamiltonian should give the centres of gravity of the full 
spin multiplets. Let us consider first the bb and cc quarkonia. For the L > 
states in our model and in many others the mass of the spin singlet compo- 
nent of each multiplet coincides with the centre of gravity of the spin triplet. 
Thus the two interpretations are equivalent. The difference occurs only for 
the S'-states. Since the success of the nonrelativistic models depends on a 
cancellation of errors, which is not understood, we cannot give a rigorous 
argument to justify one interpretation rather than the other. Let us quote, 
however, two plausibility arguments. Roncaglia and collaborators |J explain 
that the spectrum of triplet states should be particularly regular. We chose 
the centres of gravity of the triplets for the practical reasons that the masses 
of the spin singlets for the bb quarkonia are not known. A posteriori, however, 
we see that this choice has worked, which supports the conjecture that it is 
an acceptable one. For the be quarkonia the situation is even more confused, 
because spin, in general, is not a good quantum number. We assume tenta- 
tively that the eigenvalues of the nonrelativistic Hamiltonian correspond to 
the centres of gravity of the spin triplets, when the singlet-triplet mixing is 
switched off. For L > the difference between these centres of gravity and 
the centres of gravity of the whole spin multiplets are small - 3 MeV, 2 MeV 
and 0.8 MeV for the IP, 2P and ID multiplets respectively. For the S'-states 
the situation should be similar to that for the other quarkonia. 



3 



The results of our calculation for the ten bb masses of interest are given 
in Table 1. In the five cases, where comparison with experiment is possible, 
agreement is very good. A comparison with some models is also shown in 
the Table. Eichten and Quigg J7J use the potential of Buchmuller and Tye 
||. This potential is particularly respectable, because at short distances it 
reproduces correctly the two-loop result from QCD. Their agreement with 
experiment is very good for the low masses, but for the higher masses it 
deteriorates, with the error reaching almost 30 MeV for the 2P state. Kwong 
and Rosner || use the masses of the known S and P states as input, calculate 
from this input the potential and from this potential the masses of the other 
states. Thus their method avoids the bias due to a preconceived form of 
the potential. Their results for the centres of gravity of the triplets agree in 
general with ours, the greatest difference being 6 MeV for the 3P states. 

The results for the cc quarkonia are shown in Table 2. It is seen that the 
approach of Eichten and Quigg is doing much better here, about as well as 
ours. The results of Gupta and Johnson are the best. 

The results for the be quarkonia are given in Table 3. Since there are no 
experimental data, we included more theoretical predictions. Many others 
can be found from the references given in the papers used in our tables. Of 
particular interest is the comparison of our results with those of Gupta and 
Johnson [|Kj. Their fit is only slightly worse than ours for the bb quarkonia 
and slightly better than ours for the cc quarkonia. On the whole the quality 
of the two fits is good (errors not exceeding a few MeV) and comparable. 
Nevertheless, the physical assumptions behind them are very different. In 
particular Gupta and Johnson use the complete relativistic expression for 
the kinetic energy of the quarks. They also use many more free parameters 
than we do. For the be quarkonia, as shown in Table 3, the predictions of the 
two approaches differ significantly. The masses of the 3 S states calculated 
by Gupta and Johnson are larger than ours by about 40 MeV. For the IP 
state the difference is in the same direction, but smaller - about 16 MeV. 
The potential in the model of Yu-Qi Chen and Yu-Ping Kuang fJTTJ is a mod- 
ification of the potential of Buchmuller and Tye. We quote only one of the 
several versions, which they propose. The model of Roncaglia et al. |J does 
not use a potential, but obtains the masses by assuming that for a given 
kind of triplet resonances (e.g. for ground states) the mass of the particle is 
a simple function of the reduced mass of its constituent quarks. Then the 
masses of the unknown particles are obtained by interpolation or extrapola- 
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tion from the known masses. Gershtein et al. [0 use the Martin potential 
13|j supplemented with some relativistic and QCD inspired corrections. The 



scatter of the predictions is of some tens of MeV and the two models, which 
agree particularly well with the data for the bb and cc quarkonia, i.e. that of 
Gupta and Johnson and ours, are not close to each other in their predictions 
here. A comparison with the experimental data, when they come, will be, 
therefore, of great interest. 



3 Hyperfine splittings 

The spin dependent correction to the nonrelativistic Hamiltonian, which is 
responsible for the hyperfine splitting of the mass levels, is generally used in 
the form (cf. e.g. jnj) 

327T« S 1 

adapted from the Breit-Fermi Hamiltonian. The number ~ subtracted from 
the product of the spins corresponds to our assumption that the unperturbed 
nonrelativistic Hamiltonian gives the energy of the triplet. Since for the 
states with orbital angular momentum L > the wave function vanishes at 
the origin, the shift affects only the S states. Thus, the only first order effect 
of the perturbation (§) is to shift the 1 So states down in energy by 

A^ F = -^^(0)| 2 . (5) 

In order to apply this formula one needs the value of the wave function at 
the origin - this is obtained by solving the Schrodinger equation with the 
nonrelativistic Hamiltonian - and the coupling constant a s . 

Like most authors (cf. e.g. 0), we determine the coupling constant 
Oi s {m 2 c ) from the well measured hyperfine splitting of the lS(cc) state. The 
experimental value |1] 117 ± 2 MeV yields 

a s (m 2 c ) = 0.3376. (6) 

Actually, the experimental uncertainty of the measured hyperfine splitting 
introduces an uncertainty in this value, but we think that other uncertainties 
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in our calculation are more serious and we do not keep track of this particular 
uncertainty. Knowing the coupling at the scale m 2 we obtain the couplings 
at other scales as follows. The formula including the NNLO terms from [|J is 
used to correlate a s (fx 2 ) with the parameters A^ n f\ The number of flavours 
(n/) is put equal to three for jj? < m 2 (we are not interested in the region 
ji 2 < m 2 ), equal to four for > jj? > m 2 and equal to five for ft 2 > m 2 . (we 
are not interested in the region /x 2 > m 2 ). Then the value of a s {m 2 ) from 
@ is used to calculate A 1 - 3 ) and A^. Using the known value of A^ 4 ) and the 
formula from ref |IJ we find the value 

a s (m 2 b ) = 0.2064 (7) 

From this the value of A 1 - 5 ) is found and further a s (m|) is calculated. The 
value a s (m 2 z ) = 0.115 obtained from this calculation agrees very well with 
the other determinations of this parameter compiled by the Particle Data 
Group I]. Note that this supports our model, since a different choice of the 
Hamiltonian would in general lead to a different value of the wave function at 
the origin and to a different determination of a s (m 2 ) from the same hyperfine 
splitting. Then the estimate of a s {m 2 z ) would, of course, be also different. 
For the hyperfine splitting of the be quarkonia we use the coupling constant 

a s (4fi 2 bc ) = 0.2742, (8) 

so that in each case the scale is twice the reduced mass of the quark- ant iquark 
system. 

The calculated hyperfine splittings are given in Tables 1-3. No confirmed 
experimental data to check these predictions are available as yet. Let us 
note, however, that the unconfirmed experimental splitting of the 2S(cc) 
level - 92 MeV - is much bigger than expected from the models. In all cases, 
where comparison with the other models is possible, the hyperfine splittings 
predicted from our model are significantly smaller than the splittings found 
by Eichten and Quigg |7j and similar to, but usually a little larger than, the 
splittings calculated by Gupta and Johnson pTJfl . 

One can also try to compare our results with more ambitious approaches. 



A careful analysis in the framework of QCD sum rules |1J| finds for the hy- 
perfine splitting of the lS(bb) state GS-li MeV. The central value agrees very 
well with our expectation, but the uncertainty is too large to distinguish be- 



tween the potential models. A lattice calculation pi gives for the hyperfine 
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splitting of the IS (be) state 60 MeV with a large uncertainty. Again the 
central value is close to our model, but the uncertainty is big enough to be 
consistent with all the potential models quoted here. 

Let us conclude this section with two comments. The operator Hhf, 
besides shifting the 1 So energy levels by its diagonal matrix elements, mixes 
the 1 S'o states corresponding to various principal quantum numbers. Formally 
this gives significant corrections to the energy levels. These corrections, how- 
ever, are second and higher order in the perturbation. We follow the usage 
of assuming that they cancel with other second and higher order corrections 
beyond our control. At first order there are only small admixtures of other 
1 So states in any given 1 So state. These, however, seem of little interest. 
The QCD corrections to hyperfine splittings have been calculated in various 
approximations (cf. ref. jlTJ] and references given there). Since these correc- 
tions are small and controversial (it has been argued that they cancel with 
other correction [17]), we do not include them in our model. 



4 Leptonic decay widths 

The leading terms in the leptonic decay widths of the heavy quarkonia are 
proportional to the squares of the wave functions at the origin. Therefore, 
they are significant only for the 5* states. For the bb quarkonia and the cc 
quarkonia we shall consider the decays of the 3 S (vector) states into pairs of 
charge conjugated charged leptons, e.g. for definiteness into e + e~ pairs. For 
the be quarkonia we consider the decays of the 1 S (pseudoscalar) states into 
tv t pairs. Since the probability of such decays contains as a factor the square 
of the lepton mass, the decays into lighter leptons are much less probable. 

The decay widths of the vector bb and cc quarkonia into charged lepton 
pairs are usually calculated from the QCD corrected Van Royen - Weisskopf 
formula @, ]19] 



r v ^ = W4*^i ^j. (9) 

For vector mesons containing light quarks this formula leads to paradoxes 
(cf. J2(| and references contained there). For quarkonia, however, the main 



problem seems to be the QCD correction. Using the coupling constants 
a s (mg) found in the preceding section one finds that the correction linear 
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in a s is 57 per cent for cc and 35 per cent for bb. Thus in order to get 
quantitative predictions it is necessary to include higher order corrections, 
which, however, are not known. In order to guestimate the missing terms we 
tried two simple Ansatze. Exponentialization of the first correction yields 

Cfainfy) = exp ( ^-^J , (10) 

while Padeization gives 

C*Mm%)) = ' y (11) 
1 -I- — 

We use the arithmetic average of these two estimates as our estimate of the 
QCD correction factor extended to higher orders. The difference between 
C\ and C*2 is our crude evaluation of the uncertainty of this estimate. The 
resulting leptonic widths are collected in Table 4. Combining in quadrature 
the experimental errors with our estimates of the theoretical uncertainties 
we get a good overall agreement {y 2 /ND = 5.9/5). About half of the y 2 , 
however, comes from the decay width of the 2 3 S(bb), where the predicted 
value is significantly larger than the newly included experimental value M. 
Thus here there may be a problem. Let us note the relation 

9 4m^ a 2 e 2 Q 

rV->H = o T72 7 — 2\C av AE HF} (12) 

8 My a s (mQ) 

where C av is the QCD correction factor. With our choice of parameters this 
formula reduces to 

Atji 2 

IV_j, = F(Q)^AE HF , (13) 

with F(c) = 4.73 ■ 10~ 5 and F(b) = 2.33 ■ 10~ 5 . 

The formula for the leptonic widths of the pseudoscalar be quarkonia reads 

r™ = g/ijW„m;(l--g-) 2 , (M) 

where G is the Fermi constant, V c b ~ 0.04 is the element of the Cabibbo- 
Kobayashi-Masakawa matrix and the decay constant fs c is given by the for- 



mula (cf. e.g. BT 
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fk = IM^I, (15) 

Be 

where C(a s ) is a QCD correction factor. Formally this decay constant is 
defined in terms of the element of the axial weak current 

(0\A,(0)\B c (q)) = tf Bc Vc b q,. (16) 
Thus it corresponds to f w ~ 131 MeV. The QCD correction factor is fl2]J 



C{a s ) = 1 - 



be) 



7T 



m b - m c m b 
: log — 

nib + m c m c 



(17) 



With our parameters C(a s ) ~ 0.885 and since this is rather close to unity, 
we use it without trying to estimate the higher order terms. 

Substituting the numbers one finds the decay widths given in Table 4. 
The corresponding decay constants for the ground state and for the first 
excited S'-state of the be quarkonium are fg c = 435MeV and fg c = 315 
MeV. 

Let us note the convenient relation 

f 2 B = 27^ b + m 2 

which for our values of the parameters yields 



/6199 



f Bc = 57.6J—^AE HF , (19) 



where all the parameters are in suitable powers of MeV. 



5 Fine structure of the levels 

The spin dependent correction to the nonrelativistic Hamiltonian, which is 
responsible for the fine splittings, is also modelled on the Breit-Fermi Hamil- 
tonian 0, pf. It can be decomposed into a part, which is antisymmetric 



with respect to the spins of the constituents 

i v 1 / 1 1 \ / ldV(r) 8a s \ T , 
Va{t) = - -5---5T + L.(sq-Sq), (20) 
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where V(r) is the nonrelativistic potential (|l|), and a part symmetric in these 
spins. The symmetric part can be decomposed into a spin-orbit coupling 




r dr 3r 3 



1 dV(r) 8a s 



) 



L • S + 



3m Q mQ r 



4a 



s 



U-S, (21) 



where 



S = S Q + S Q' 



(22) 



and a tensor part 




4a s 1 
3 r 3 



[3(sq -n)(sQ -n) -s Q -Sq] 



(23) 



where the versor n = -. In the first perturbative approximation these cor- 
rections to the Hamiltonian not only shift the mass levels, but also mix some 
states with different values of the orbital angular momentum and spin. 

Let us note some useful selection rules. Angular momentum and parity 
are good quantum numbers, therefore states with different J and/or P do not 
mix. States with orbital angular momenta differing by one unit do not mix, 
because parity is a good quantum number, and states with angular momenta 
differing by more than two units cannot mix because of the Eckart-Wigner 
theorem. From symmetry with respect to the exchange of spins in the LS 
basis, the operator Va can only contribute to matrix elements between the 
spin singlet and spin triplet states. The other two operators contribute only 
to the matrix elements between spin triplet states. When Q is the charge 
conjugate of Q, charge conjugation C = (— 1) L+5 is a good quantum number 
and states with different spins do not mix. The matrix elements between 
spin singlet states vanish; thus, when spin is a good quantum number, we 
predict within each spin multiplet with L > : 



Here the centre of gravity of the triplet is defined by the usual formula 



MCL)=M c . . g fL) 



(24) 



M C . . 9 .( 3 L) 



(2L + 3)M( 



3 L i+1 ) + (2L + 1)M( 3 L L ) + (2L - 1)M( 3 L L ^) 



3(2L + 1) 



(25) 
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There is no firmly established experimental data to test this prediction, but 
most models and the preliminary data for the l 1 Pi(cc) state agree with it 
within about f MeV. Of course, this prediction is common to all the models, 
which use the spin dependent corrections to the Hamiltonian as given here. 

Since the wave function of the quarkonium in the LS basis factorizes into 
a space part and a spin part, the matrix elements of the space operators 
and of the spin operators can be calculated separately. The necessary matrix 
elements for the spin operators between spin triplet states are 



(J, L', 1|L • S| J, L, 1) = i[J(J + 1) - L(L + 1) - 2}5 LU (26) 

(L + 1, L, l|3(s Q • n)(sg • n) - s Q ■ Sq\L + 1, L, 1) = 

2L 



2L + 3' (27) 
(L, L, l|3(s g • n)(sg • n) — Sq • Sq|L, L, 1) = 2, (28) 

(L - 1, L, l|3(s Q • n)(sg • n) - s Q • s^\L - 1, L, 1) = 



2(L + 1) 
2L-1 ' 

(L + 1, L - 2, 1|3(sq • n)( S Q • n) - s Q • Sq|L + 1, L, 1) = 



3 JL(L - 1) 



(29) 



(30) 



2(2L - 1) 

The only necessary nonzero matrix element between a spin singlet and spin 
triplet state is 



(L, L, 0|L • (s Q - Sq)|L, L, 1) = y L(L + 1) (31) 

Let us consider now the two matrix elements in coordinate space. The cal- 
culation of the matrix elements of the operator l^YSll [ s standard, but the 
calculation of the matrix elements of the operator % requires an interpreta- 
tion of a s . We propose to interpret a s as a function <5 s (r) defined as follows 

~ / n 12 tt (A (n 'V) 2 - 1 

= oo o r UTT ^ , o1 > ( 32 ) 



33 - 2n/ log [(A( n /)r) 

where rif equals three for r < equals four for ^- < r < ^- and equals 
five for r > The parameter A( 4 ) is obtained from the conditions a(-^-) = 
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a(m^) and <5(^-) = a(m^). Each of these conditions gives a slightly different 

A^ 4 ). We use the geometric mean of the two results. They are so close to 
each other that taking the arithmetic mean instead of the geometrical one 
makes no difference within our precision. Knowing we recalculate a s (r) 
at r = — and r = — and fix A^ 3 ^ and A^ 5 ^ so that the function aJr) is 
continuous at these points. We find 

A (3) = 0.f657 GeV, A {4) = 0.1384 GeV, A (5) = 0.1015 GeV. (33) 

Our form of the function a s (r) is, of course, inspired by the standard one- 
loop formula for a s . The numerator is introduced in order to compensate 
the zero of the denominator at A^ 3 V = 1. Its exact form has little effect in 
the range of r dominating the integrals. 

The calculations for the bb and cc quarkonia, where C is a good quantum 
number, and for the J = L ± 1 states of the be quarkonia, which must be 
spin triplets, involve only the energy shifts due to the symmetric spin orbit 
interaction and to the tensor interaction. The singlet and triplet J = L states 
of the be quarkonia mix under the influence of the antisymmetric spin-orbit 
interaction. The results of the calculations are given in Tables 1-3. For the bb 
quarkonia experimental data is available for the splittings of the IP and 2P 
states. Agreement between our model and this data is within about 1 MeV. 
The agreement with other models is within 5 MeV, i.e. much better than 
for the hyperfine splittings. A similar agreement with the model of Kwong 
and Rosner holds for the 3P states, but for higher angular momenta the 
discrepancies increase. We predict much larger splittings. In particular for 
the F states we predict a splitting of about 13 MeV, while Kwong and Rosner 
expect a negligible splitting within about 1 MeV. Even for the G states we 
expect a splitting of about 10 MeV. Thus, the L-dependence of the fine 
splittings is seen as an important observable to distinguish between models. 
For the lP(cc) states our splittings agree with experiment and with the very 
good predictions of Gupta and Johnson about as well as for the centres of 
gravity of the triplets i.e. within about 5 MeV. The splittings predicted by 
the model of Eichten and Quigg are too small by about a factor of two. 
For the be quarkonia our predictions for the IP states agree with Gupta 
and Johnson within about 2 MeV except for J = 0, where our splitting is 
smaller by 10 MeV. For the ID and 2P states there are only the predictions 
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of Eichten and Quigg and of Gershtein et al. [12] for comparison. There 
is rough qualitative agreement except for the 3 Pi state, where we predict a 
down shift by almost 20 MeV, while the other models find only very small 
shifts. 

The mixing between the spin singlet and the spin triplet states can be 
parameterized in terms of mixing angles. We find sin^ip = 0.374, sin02P = 
0.385 and sin0i£> = 0.244. Thus the mixing within the two P multiplets is 
almost the same, while the mixing among the D states is somewhat smaller. 



Both these results contradict the results of Gershtein and collaborators ]|i2 
who find that mixing increases when going from the IP to the 2P states 
and from the 2P to the ID states. A possible reason for this discrepancy is 
that these authors use for the mixing formulae, which are different from ours. 
In particular their singlet-triplet mixing does not vanish for tuq = m-Q. As 
compared with Eichten and Quigg, who have calculated mixing only for the P 
states, we have rough agreement for the 2P states (they find sin 2 p = 0.290), 
while they find much less mixing for the IP states (sin^ip « 0.06). 

The mixing of spin triplet states differing by two units of orbital angular 
momentum (L — 2, L) is small. In our model we find mixing angles of order 
10~ 3 or less. It seems of little interest, except that it enhances the leptonic 
decay widths of the L > 2 states (cf. e.g. [p2 ). This enhancement, however, 



is difficult to calculate reliably, because a given high L state mixes with 
various L — 2 states and the states above the strong decay threshold are also 
important for this analysis. 



6 Conclusions 

We propose a model containing six free parameters: the three parameters in 
the nonrelativistic potential ([[]), the masses of the c and b quarks @ and 
the strong coupling at the m c scale (|6|). This model is applicable to all the 
heavy quarkonia below their strong decay thresholds. 

We obtain for the hh quarkonia 12 quantities (five spin averaged masses, 
four independent mass differences due to fine splittings and three leptonic 
decay widths) in good overall agreement with experiment. The least suc- 
cessful predictions are for the fine structure shift of the 2 3 P X state, which 
is measured to be —4.8 ± 0.5 MeV, while we find —6.0 MeV and for the 
leptonic decay width of the 2 3 S(bb) state, where we find (0.59 ± 0.03) KeV, 
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while the newly included experimental result ffl] is (0.52 ±0.03) KeV. For the 
cc quarkonia we find 8 quantities, which can be compared with experiment 
(six masses and two leptonic widths). Here in most cases the difference be- 
tween the measured value and the prediction exceeds the experimental error, 
but the errors in the mass predictions do not exceed a few MeV and the 
errors in the leptonic decay widths do not exceed 1.6 s.d.. On the whole, 
with 6 parameters we predict 20 quantities in good (for bb) or fair (for cc) 
agreement with experiment. The only other model known to us, which has a 



comparable record, is the model of Gupta and Johnson fL0| , but this model 
has many more free parameters. 

We give predictions for the yet unmeasured masses of the quarkonia and 
for the leptonic widths of the be quarkonia. The predictions are listed in 
Tables 1-4. Here we would like to make the following general remarks. Our 
model predicts much larger fine splittings at high L than the model of Kwong 
and Rosner ||. We also find significantly heavier be resonances than Gupta 
and Johnson [TTJ, which is remarkable, because the two models give similar 



descriptions of the charmonia and of the bottomonia. The discrepancy is in 
the spin averaged masses. The splittings of the levels, except for the hyperfine 
splitting of the IS* level, are similar. 
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Table captions 

Table 1. Mass spectrum of the bb quarkonia below the threshold for strong 
decays (2tob = 10558 MeV) in MeV. AX denotes the difference between the 
mass of particle X and the centre of gravity of the spin triplet part of the 
multiplet, where X belongs. 

Table 2. Mass spectrum of the cc quarkonia below the threshold for strong 
decays (2m/) = 3729 MeV) in MeV. AX denotes the difference between the 
mass of particle X and the centre of gravity of the spin triplet part of the 
multiplet, where X belongs. 

Table 3. Mass spectrum of the cb (be) quarkonia below the threshold for 
strong decays (m D + m B = 7143 MeV, m B * + m D = 7189 MeV) in MeV. 
AX denotes the difference between the mass of particle X and the centre 
of gravity of the spin triplet part of the multiplet, where X belongs, (a) - 
particle above its strong decay threshold. 

Table 4. Leptonic widths in keV. 
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